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ABSTRACT 

We present 12 CO J = 1 — > observations of ultraluminous infrared galaxies (ULIRGs) obtained 
using the Redshift Search Receiver (RSR) on the 14-m telescope of the Five College Radio Astronomy 
Observatory. The RSR is a novel, dual-beam, dual-polarization receiver equipped with an ultra- 
wideband spectrometer backend that is being built as a facility receiver for the Large Millimeter 
Telescope. Our sample consists of 29 ULIRGs in the redshift range of 0.04-0.11, including 10 objects 
with no prior 12 CO measurements. We have detected 27 systems (a detection rate of 93%), including 
9 ULIRGs that are detected in CO for the first time. Our study has increased the number of local 
ULIRGs with CO measurements by ~ 15%. The CO line luminosity L' co , correlates well with far- 
infrared luminosity Lfir, following the general trend of other local ULIRGs. However, compared to 
previous surveys we probe deeper into the low CO luminosity end of ULIRG population as a single 
study by including a number of CO faint objects in the sample. As a result, we find 1) a smoother 
transition between the ULIRG population and local QSOs in Lfir — L' co ( "star formation efficiency" ) 
space, and 2) a broader range of Lfir/L' co flux ratio (~60-10 3 L Q /[K km s _1 pc 2 ]) than previously 
reported. In our new survey, we also have found a small number of ULIRGs with extreme Lpm/ L' co 
which had been known to be rare. The mid-/i? color and radio-excess of 56 local ULIRGs as a function 
of F I R-to-CO flux ratio is examined and compared with those of spirals/star burst galaxies and law-z 
QSOs. In this paper, using a large sample of local ULIRGs we explore the origin of their current 
power source and potential evolution to QSOs. 

Subject headings: galaxies: evolution — galaxies: starburst — ISM: molecules 



1. INTRODUCTION 

An important discovery from the IRAS all-sky sur- 
vey is a large population of infrared luminous (Ljr > 
10 11 Lq) galaxies that are emitting the bulk of their lu- 
minosity in the infrared (Sanders & Mirabel 1996). The 
most extreme objects, so-called "ultraluminous infrared 
galaxies (ULIRGs)" with Ljr > 10 12 Lq, are the most 
luminous objects in the local Universe (Solomon et al. 
1997). ULIRGs are rare at low redshift (z < 0.3), but 
their co-moving IR energy density grows rapidly with 
increasing redshift (LeFloc'h et al. 2005), and galaxies 
with ULIRG-likc luminosity are quite numerous at z > 1 
(Farrah et al. 2003, and references therein). ULIRGs 
therefore represent an important population for under- 
standing the cosmic star formation history and the cos- 
mic energy budget. 

Most ULIRGs appear to be advanced mergers (Sanders 
et al. 1988a; Farrah et al. 2001) while -15% of IR lu- 
minous systems show characteristics of quasi-stellar ob- 
jects (QSOs) in their optical spectra (Sanders & Mirabel 
1996). Whether the origin of their large IR luminosity 
is a merger-driven starburst or AGN activity (or both) 
is still not well established. Their eventual fate is also 
controversial. In the multi-wavelength study of a flux- 
limited sample, Sanders et al. (1988a, b) have suggested 
that ULIRGs hosting an AGN are transition objects that 
eventually evolve into optically selected QSOs as the 
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AGN light dominates the decaying starburst light. 

A common trait found among the ULIRG population 
is the large molecular gas content (M# 2 > 10 10 M©), in- 
ferred from their observed CO luminosity (Sanders et al. 
1991; Young et al. 1995; Solomon et al. 1997), capable of 
fueling both the bursts of star formation and the AGN 
activities seen in these systems. Indeed, the correlation 
between the far-infrared luminosity Lfjr and the CO 
line luminosity L' co of ULIRGs forms a continuous track 
with those of QSOs in the Local Universe (Riechers et al. 
2006, and references therein), which is supportive of the 
proposed ULIRG-QSO evolution scenario. QSOs with 
LiR > 10 12 Lq however are rare in the Local Universe, 
implying that ULIRGs are not simply a stage of transi- 
tion to QSOs (Farrah et al. 2001; Yun et al. 2004) or that 
the transition may occur quite rapidly (Yun et al. 2004). 

Considering the rich history of multi- wavelength inves- 
tigations of the ULIRG phenomenon and their impor- 
tance as the local analogs to the luminous dusty galax- 
ies found in the early universe (so-called "submillimc- 
ter galaxies (SMGs)"; see review by Blain et al. 2002), 
it is surprising that the total number of published CO 
measurements for ULIRGs is limited to only —50 or so, 
primarily from the surveys conducted by Sanders et al. 
(1991) and Solomon et al. (1997) 3 , and by Mirabel et al. 
(1990) for southern objects. There have been also some 
CO studies of individual ULIRGs in the Local Universe 
(z < 0.3, e.g. Combes et al. 2006), although these do not 
contribute to the total number statistics by large. 

3 More recent ULIRG surveys of the higher J CO transitions 
(Yao et al. 2003; Narayan et al. 2005) have targeted mostly the 
same sample. 
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This limited information on CO emission and molecu- 
lar gas content reflects both the rarity of ULIRGs in the 
local volume and the challenge of detecting CO emis- 
sion from galaxies even at these modest distances. It 
raises a real concern that the existing CO data may be 
too small to characterize the molecular gas properties 
of the ULIRG population with sufficient statistics. The 
existing CO data may also reflect strong biases in the 
sample selection and/or the survey strategies adopted 
by these earlier studies. In comparison, the total num- 
ber of CO measurements available for the z > 1 SMGs 
is >30 systems (see review by Solomon & vanden Bout 
2005), which is generally considered too few to draw any 
broad conclusions about the molecular gas content and 
properties of the SMC population. 

In order to address these issues, we have conducted 
a new CO survey of 29 ULIRGs in the redshift range 
z « 0.04-0.11 using the sensitive, ultra-wideband spec- 
trometer system called the Redshift Search Receiver 
(RSR) on the Five College Radio Astronomy Observa- 
tory (FCR AO) 14-m telescope. The broad spectral cover- 
age of the RSR enables a CO line detection even in cases 
where the reported optical redshifts might have been sig- 
nificantly in error. Our sample is comparable in size to 
the largest previous CO surveys of ULIRG population 
(e.g. Solomon et al. 1997), but probes deeper into the 
CO luminosity by including more CO faint objects in 
the sample as a single study. Our sensitivity allowed us 
to detect a number of sources with the comparable line 
intensity as the faintest ULIRGs in CO that have been 
known to date (e.g. IRAS 08572+3915 by Solomon et al. 
1997 or IRAS 10173+0828 by Sanders et al. 1991). 

This paper is organized in the following order. In §2 we 
describe our sample selection criteria and the properties 
of the sample. In §3 we describe the observations and 
data reduction. In §4 we present the results and discuss 
the data quality. In §5 the origin of IR emission and 
the evolution of the local ULIRGs are discussed. Finally 
in §6 we summarize our findings. We assume a stan- 
dard ACDM cosmology with H =71 km s _1 Mpc -1 , 
Q M =0.27, and Q A =0.73 (Spergel et al. 2003) through- 
out the paper. 

2. SAMPLE 

Our sample is constructed primarily from the 1 Jy 
ULIRG 4 sample, which consists of 118 sources identi- 
fied from the IRAS Faint Source Catalog (Moshir et al. 
1990) with Seon > 1 Jy, 5 > -40°, and \b\ > 30° (Kim 
& Sanders 1998). All 18 ULIRGs in the 1 Jy Sample 
whose 12 CO J = 1 ^ line falls within the frequency 
range of J^bs « 104-111 GHz (z =0.043-0.11, see §3) were 
selected as the primary sample. In addition, we have se- 
lected four targets from the 2 Jy sample of Murphy et 
al. (2001) with S 60f > > 1.94 Jy, <5 > -35°, and |6| > 5°. 
Lastly, in order to provide comparison with our obser- 
vations (§ 4.2), 20 ULIRGs with previous CO measure- 
ments were added from Mirabel et al. (1990), Sanders et 
al. (1991), Solomon et al. (1997), and Gao & Solomon 
(1999). Those ULIRGs were all previously detected ex- 
cept for IRAS 20414-1651 (Mirabel et al. 1990). These 

4 The conventional definition of an ultraluminous infrared galaxy 
is a galaxy with rest-frame 1-1000 £tm luminosity of Lm > 10 12 Lq 
(see Sanders & Mirabel 1996). 



include 13 sources which overlap with the sample from 
Kim & Sanders (1998) or Murphy et al. (2001), or both, 
yielding the total number of our sample 29. General 
properties of the 29 ULIRGs are summarized in Table 1. 

The optical Digitized Sky Survey plate images of all 
29 ULIRGs and the color images of 12 available ULIRGs 
from the Sloan Digital Sky Survey data archive are 
presented in Figure 1 and 2. As shown in the opti- 
cal images, morphological peculiarities are quite com- 
mon in ULIRGs. Many galaxies show faint stellar tails 
(e.g. IRAS 12112+0305), rings (IRAS 15250+3609), or 
extended/ asymmetric low surface brightness disks (e.g. 
IRAS 10565+2448). More than 75% of the sample are 
classified as Hii, LINER or a starburst systems. Three 
objects IRAS 05189-2524, IRAS 08572+3915, and 
IRAS 12540+5708 (Mrk 231) have previously been iden- 
tified as "warm" ULIRGs with /2s //60 > 0.2 by Sanders 
et al. (1988b), and another three (IRAS 14394+5332, 
IRAS 15130-1958, and IRAS 15250+3609) have warm 
IR color (0.18 < /25//6O £ 0.2), which is indicative of 
active galactic nuclei (AGN; de Grijp et al. 1985). In 
addition, IRAS 00057+4021 and IRAS 15462-0450 are 
classified as Seyfert type 2 and 1, respectively, and likely 
harbor an AGN. 

Our study is one of the most systemic and deepest 
CO surveys of a fairly large sample of ULIRGs to date, 
including 9 objects with no prior CO measurments and 
one which was observed but not detected. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. Observations with the RSR on the FCR AO H~m 

Telescope 

The observations presented here were conducted us- 
ing the Five College Radio Astronomy Observatory 
(FCRAO) 14-m telescope between March and May in 
2007 and in 2008 May. The FCRAO 14-m telescope is 
a radome-enclosed single-dish millimeter telescope, lo- 
cated on Prescott peninsula within the Quabbin Reser- 
voir, Massachusetts. The spectra were taken using the 
Redshift Search Receiver (RSR). The RSR is a sensitive, 
ultra-wideband spectrometer that is being built at the 
University of Massachusetts as one of the facility instru- 
ments for the 50-m diameter Large Millimeter Telescope 
(LMT). 

The RSR consists of a frontend receiver which uses 
a novel construction for a millimeter wavelength system 
and a set of wide-band analog auto-correlation spectrom- 
eters. There are four receivers each covering 74-111 GHz 
instantaneously in a dual-beam, dual-polarized system. 
The input includes a novel electrical beam switch which 
operates at 1 kilohertz (kHz) to overcome the 1/f noise 
originating within the frontend amplifiers as well as at- 
mospheric noise to ensure excellent baseline stability. 
The frontend uses MMIC (monolithic microwave inte- 
grated circuit) amplifiers and two very wideband mixers 
to convert each receiver band to two intermediate fre- 
quency (IF) channels. After further conversion, the IF 
signal passes into a spectrometer based on analog auto- 
correlation. Sets of tapped delay lines sample and multi- 
ply the signal with progressive delays to generate a spec- 
trum with 31 MHz resolution. Six spectrometers, each 
covering 6.5 GHz bandwidth, are used with each pixel. 
The entire 36 GHz of bandwidth of each pixel is handled 
by six boards, each having 256 lags (in all 1536 lags per 
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TABLE 1 

General Properties of the Sample of 29 ULIRGs a 
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a These data have been compiled from the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National Aeronautics and Space Administration. 

b References for previous CO measurements: M90-Mirabel et al. (1990); S91-Sanders et al. (1991); S97-Solomon et al. (1997); GS99-Gao & 
Solomon (1999) [References for the optical studies: KS98-Kim & Sanders (1998)'s Uy sample; MOl-Murphy et al. (2001)'s 2 Jy sample] 



pixel). Occasionally, a few of the lags develop problems, 
which arc flagged and blanked out before converting to 
the frequency domain. However, less than 1% of all lags 
had such problems, and most of these faulty lags have 
since then been fixed in the spectrometer when the re- 
ceiver was taken off the telescope More details of the RSR 
and the LMT can be found in Erickson et al. (2007). 

During the observing season in 2007, all four frontend 
pixels (receivers) were available while only 4 of the 24 
backend spectrometers had been fabricated. This yielded 
a set of 4-frontend pixels x 1-backend spectrometer with 
a bandwidth range of 104-111 GHz (1st LO could be 
tuned) , and this spectral coverage ultimately limited the 
rcdshift range in our sample selection. In 2008, 12 spec- 
trometers became available, and we had a set of 4-pixels 
x 3-spectrometers with an instantaneous frequency cov- 
erage of 92-111 GHz. Therefore some of our new spectra 
cover a much larger spectral (redshift) range, but we fo- 
cus on the analysis of only the 12 CO J = 1 — > transition 
in this paper. 

We integrated between 1 and 11 hours (Table 1) on 
each source depending on the weather conditions and 
the peak of the CO intensity. The system temperature 
during the run varied from 165 to 380 K with a typical 



value of T sys ~ 200 K in both years. For calibration, 
normally we obtained one 5-min scan on the reference 
sky every three consecutive 5-min scans on the target. 
When the weather was bad (T sys > 400 K), we did cal- 
ibration every other scan although most scans obtained 
under those weather conditions were thrown away in the 
final co- addition. As a result, we achieved a typical rms 
sensitivity of a < 0.5-0.6 mK in T\ (corrected antenna 
temperature for the atmosphere absorption and spillover 
losses). We checked the pointing and the focus using 
planets every 2 to 3 hours. The antenna gain was moni- 
tored by measuring the intensity of planets. An average 
gain factor of G= 45 Jy K _1 is adopted (in T\ scale). 



3.2. Data Reduction Using SPAPY 

The RSR data have been reduced using the SPAPY (SPA 
+ PYthon). SPAPY is a data reduction software that has 
been developed by G. Narayanan mainly for the RSR 
data reduction. For each galaxy, the final spectrum is 
constructed using the following procedure. First a set 
of scans of different pixels were individually inspected in 
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Fig. 1— Optical images from the 2nd Digitized Sky Survey (DSS2) of the sample of 29 local (z = 0.043-0.11) ULIRGs. Each field 
(1.5' X 1.5') is centered on (a, 5)2000 as shown in Table 1. The IRAS catalog number is indicated on the top of each panel. A circle of 50" 
(sibeamsizc of 14m at 3mm) is shown in solid line. 



order to exclude data obtained under bad weather con- 
ditions or with serious instrumental noise. Then con- 
tinuum channels are selected that exclude channels with 
CO line emission and the end channels affected by the 
bandpass roll-off. A linear baseline is fitted over these 
continuum channels to calculate the rms noise of each 
scan. A set of scans for each object are summed weighted 
by the rms noise. Lastly a baseline is removed from 
the combined spectrum, using the same frequency range 
where we measure the rms of individual scans. For most 
galaxies, a linear baseline has been fitted while a second 
order fit has been adopted in some cases, especially for 



galaxies with low signal-to-noisc. The second order fit 
barely changes the rms over the entire band (< 5%) but 
it helps to bring out weak line features by modifying the 
local baseline structure around the CO line. The mea- 
sured rms value of the fully processed spectrum for each 
galaxy is listed in Table 1. 

In Figure 3 (a)-(d), the coadded spectra of 29 ULIRGs 
are shown, which have been divided into four groups de- 
pending on the signal-to-noise. On the left-side of each 
object, the full band spectrum, the full band spectrum 
excluding noisy end channels is presented in the same T* A 
scale. Twenty seven objects in group (a), (b), and (c) are 
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Fig. 2. — The optical color images from the Sloan Digital Sky Survey (SDSS, http://www.sdss.org/) of 12 ULIRGs. The size of each 
field is the same as Figure 1. 



qualified as detection by revealing a > 3tr feature within 
a few 100 km s _1 around the inferred CO frequency from 
the optical redshift, which are marked with dotted lines 
in the full band spectra in Figure 3 (a)-(d). Except for 
2 ULIRGs in group (d) which fail our definition of de- 
tection, the rest 27 have been classified from (a) to (c) 
based on the following criteria: 

(a) Bright or Intermediate: In seven ULIRGs with 
S/N> 7 (bright), and in another seven objects with 
3.5<S/N< 7, Wco > 250 km s" 1 and the peak flux 
density F]~ 1 Q k > 1.5 mK (intermediate), the lines 
can be picked out without any difficulty. 

(b) Faint or Narrow: Among eight ULIRGs with 
3.5<S/N<5.5, ve objects with F^ k < 1.5 mK but 
Wco >200 km s" 1 (faint), and three objects with 
Wco <200 km s" 1 yet F^ k > 1.4 mK (narrow) 
still show well dened line features. 

(c) Marginal: There are five ULIRGs that still meet our 
detection criteria but the lines are not as prominent 
as the other two groups (1< F^ k < 1.2 mK). There 
could be some confusion with nearby noise structure 
as indicated by the error in Wco of > 40%. 

As seen in the full band spectra, the baselines are fairly 
flat throughout the entire frequency coverage with well 
behaved noise. This is also demonstrated in Figure 4 



where the rms from the full band (Table 1) are com- 
pared with the local rms which were measured using the 
continuum channels within ~ 1 — 2 GHz range around 
the CO line. In most cases, the rms are found along the 
equivalent line with < 10% of variation. 

It should be noted that such flat baselines over such 
wide bandwidths have hitherto not been reported for 
millimeter wave receivers, and is a result of the inno- 
vative design of the RSR. The flat baseline of the RSR 
with well behaved noise will be particularly important 
for a blind line search of galaxies with no previous red- 
shift measurement, and hence this demonstration has 
been one of our goals in this study. These character- 
istics of the RSR also allowed us to integrate down to 
detect fainter lines, yielding 93% of the detection rate 
including 5 marginal detections. The two non detections 
might be due to the insufficient band coverage. In partic- 
ular, Gao & Solomon (1999) reported the CO luminosity 
of IRAS 16487+5447 which is almost three times larger 
than our upper limit based on the CO line intensity up- 
per limit (3xrmsx250 km s _1 , §4.1). However, our base- 
line fit has successfully brought up the line close to the 
band edge even in some cases with marginal signal-to- 
noise. Therefore we suggest that those two galaxies may 
contain molecular gas of low flux density with a huge line 
width (> 500 km s -1 ) as found in some extreme cases 
like IRAS F11180+1623 by Dinh-V-Trung et al. (2001). 

4. RESULTS 
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Fig. 3. — (a) The RSR spectra of 14 ULIRGs with a high (left column) to intermediate (right column) CO strength. The quantitative 
descriptions for the classication are given in §4.2. On the left side of each ob ject, the full abdn spectrum is shown except for a few 
noisy channels at the edge. On the right-hand side, each spectrum is zoomed-in around the CO emission. The full-band spectra with the 
frequency range of 103.7-110.9 GHz are shown in the same T a (antenna temperature corrected for the atmosphere absorption and spillover 
losses) scale for all galaxies. The dotted lines in each full-band spectrum represent the frequency range used to measure the CO intensity 
and derive the CO line width. The upper arrow indicates the estimated CO frequency from the optical redshift. The zoom-in spectra of 
1.2 GHz (~3100 km s _1 ) width centered on the CO frequency are scaled down with the peak of the CO emission. The dotted line in the 
zoom-in spectrum represents the centroid CO frequency. 



4.1. Measuring the CO Quantities 

Integrated intensity (/co), center frequency (yco) and 
CO line width (Wco) are determined from the final re- 
duced spectra. The integrated line intensity in K kms -1 
is calculated by summing the flux density of the line 
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Fig. 3. — (b) The RSR spectra of 5 ULIRGs that are faint in CO (left column) and 3 ULIRGs that are faint and narrow (right column). 
See the caption of Figure 3- (a) for more details. 



where Fu ne is flux density in K at each line channel, Av 
is the channel width (0.031 GHz for the RSR), vu n e is 
the central frequency of each spectral channel in GHz, 
and c is the speed of light in km s _1 . The central fre- 
quency, vqq has been calculated as a first moment of the 
spectrum in the same frequency range used to measure 
Ico- The line width, lUco, is calculated as /co/^peak- 
For the two non-detections, the quoted 3cr upper-limits 
are estimated assuming a line width of 250 km s _1 
(cf. the mean line width Wco of the 27 CO detected 
galaxies^ 264 ± 60). The rms noise locally measured 
around the CO line (< 2 GHz, Table 2) has been adopted 
to calculate the upper limits. The measured CO quanti- 
ties are summarized in Table 2. 

4.2. Comparisons with Previous Measurements 

Among the prior CO studies of LIRGs/ULIRGs, 
Solomon ct al. (1997) collected coherent CO data us- 
ing the IRAM 30m telescope of 37 ULIRGs and is the 
most comparable with our study. A comparison of our 
Ico of 1 1 ULIRGs which were also observed by Solomon 
et al. (1997) is shown in Figure 5. Gain factors G =45 
and 4.5 Jy K _1 in TjJ and T TO & scale are adopted for 



the 14m and 30m telescope, respectively. In Table 3, 
Ico from the 14m telescope is given in Jy km s -1 . As 
shown in Figure 5, most of the objects are within their 
la measurement uncertainty from the /co,30m = ^co,i4m 
line and Ico from the two telescopes are overall in good 
agreement. 

In Figure 6, we directly compare the 14m spectra with 
those from the 30m for four objects with a range of signal- 
to-noise: one CO bright object (112540, also known as 
Mrk 231), one with intermediate brightness (100057), one 
CO faint object (110494), and one with a marginal de- 
tection (108572). As shown on the left side of the figure, 
the 14m spectra are generally in good agreement in spite 
of the low spectral resolution (^90 km s _1 vs. 24 or 
48 km s^ 1 for 30m). There may be some confusion in the 
cases of marginal detections although the discrepancy in 
the total flux is expected not more than 40% inferred 
from the comparison of 108572. Lastly, the rcdshifts de- 
termined from our CO spectra (Table 3) are consistent 
with optical redshifts for the 27 detected sources, with 
A<V CO - V opt > = 121 ± 93 km s" 1 . 

A comparison of the distribution of CO luminosities 
for our survey sample and those in the same redshift 
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Fig. 3. — (c) The RSR spectra of 5 ULIRGs of marginal detections. See the caption of Figure 3-(a) and §3.2 for further details. 



Not detected 




104 105 106 107 108 109 110 

Frequency (GHz) 



103.5 104 104.5 



I ■' 1 :l I 1 1 1 I 
116487 



I 'i I. | I | I ... I ... I ... I ... I 



104 105 106 107 108 109 110 
Frequency (GHz) 



104 104.5 105 



Fig. 3. — (d) The RSR spectra of 2 ULIRGs which were not detected in our study. The dotted lines in the full-band spectra indicate the 
1 GHz width (~ 2600 km s _1 ) around the inferred CO frequency the optical velocity. The y-axis in the zoom-in spectrum corresponds to 
—3(7 to 5c as presented in Table 1. 



range measured by Solomon ct al. (1997) is shown in 
Figure 7. This plot nicely shows that nearly all of the 
ULIRGs in our sample are in the lower half of the CO 
luminosity bins. Probing deeper into the lower CO lu- 
minosity objects, we find a broader range of L fir/ Leo 
ratio for the ULIRG population, including the discovery 
of a handful of ULIRGs with extremely large Lfir/ Leo 
ratios (> 250 Lq/L{). This newly identified population 
is discussed further in § 5.1. 

4.3. CO Luminosity and Molecular Gas Mass 
The CO line luminosity, L' co , is calculated as 

L' co = 3.25xl0 7 W~ 6 2 s £>| (l+z)" 3 K km s' 1 pc 2 (2) 

where Iqo is the velocity integrated CO intensity in 
Jy km s _1 , Dl is the luminosity distance in Mpc, v a b s is 
the observed frequency in GHz, and z is the redshift of 
the object (Solomon et al. 1992). We used the redshift 
determined using the CO line to calculate Dl- For the 
undetected sources, L' co upper limits are computed us- 
ing the upper limits of Ico m Table 2 and their optical 
redshifts. Derived CO luminosities are listed in Table 3. 

The CO luminosity of our sample varies by more than 
a factor of 10, ranging 1.2 - 15.3xl0 9 K km s _1 pc 2 . 
Adopting the H 2 -to-CO conversion factor for starburst 



systems {M gas /L[ 



CO 



0.8 Mq [K km s- 1 pc 2 ]- 



Downcs & Solomon 1998), L' co of our sample corre- 
sponds to 1 - 12xl0 9 Mq of H 2 , which is comparable 
to normal to gas-rich spiral galaxies (Gao & Solomon 
2004b). The median L' co of our sample, 6 — —7 x 
10 9 K km s _1 pc 2 , is about three times as large as the L* 
value of the local CO luminosity function (Keres et al. 
2003), consistent with a simple model that these ULIRGs 
are mergers of two massive spiral galaxies. 

5. DISCUSSION 

5.1. Non-linear Correlation between Lfir and L' co 

Using the fluxes in the IRAS Faint Source Catalog, we 
have calculated FIR luminosity, Lfir following Lons- 
dale et al. (1985), 

J^FIR = 4nD 2 L FIR L Q (3) 
Dl is the luminosity distance (see § 4.3) and FIR is the 
far-infrared flux, 

FIR = 1.26 x 10~ 14 (2.58/60 + fwo) W mT 2 \ (4) 
where /go and /inn are flux densities at 60 and 100 /im 
in Jy, respectively. Considering the large uncertainty in 
dust emissivity, we have not applied any color corrections 
to L fir in Table 3. 

There is a well-known correlation between IR lu- 
minosity and CO luminosity for star forming galax- 
ies, presumably driven by the scaling relation between 
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Fig. 5. — Comparisons of Iqq from the FCRAO 14m telescope 
with those measured by the IRAM 30m telescope (Solomon et al. 
1997). The antenna gain, G=45 and 4.5 Jy K _1 in T' A and T mb 
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Star formation activity and the amount of fuel available 
(see review by Young & Scoville 1991). This correla- 
tion is reproduced in Figure 8, including the 29 new 
CO measurements of ULIRGs obtained using the RSR 
system. Galaxies with normal to high star formation 
rate (Solomon et al. 1997; Gao & Solomon 2004b) and 
Palomar-Green (PG) QSOs (Alloin et al. 1992; Evans et 
al. 2001; Scoville et al. 2003) are also shown for compar- 
ison. As previously reported (e.g. Sanders et al. 1991; 
Solomon et al. 1997; Gao & Solomon 2004b), ULIRGs 



TABLE 2 
CO J = 1 — > Line Parameters 



\J DJCCt 


rms 






Wco 




mK 


K km s — 1 


GHz 


km s" 


-i 


1UUUO / 


ft CI 
U.Ol 


n snj-n 1 o 
U.oUztU.lz 


1 1 n qk j_n no 
llU.oOztU.Uz 


347± 


91 


lUOUoo 


U.Ol 


i /i n_Ln i o 
i .4yztu. iz 


i nn qc j_n m 
iuy .ooztU.Ul 


177± 


20 


mci QO 
lUoloy 


n sq 

U.OO 


n fifij_n 1 
U.OOztU.10 


1 1 n j_n c\a 
11U.OO±U.U4 


258± 


93 


TH2 C70 

lUoO i z 


n Q/l 
U.o4 


n qi -i~n ns 
U.olztU.Uo 


i ns Qfi j_n no 
lUo.oOztU.Uz 


290±116 


lUyill 


U.OO 


n ^ e j_n 1 o 
U.OoztU.lz 


inn Q9J.n H7 
lUy.ooztU.U I 


302± 


94 


n nnQ c 
llUUoo 


n en 

u.oy 


j_n 1 q 
1.14±U.lo 


i nc QfiJ-fi m 
lUo.oUztU.Ul 


320± 


63 


T1 m 7Q 
11U1 i o 


n qq 
U.oo 


n ofi-Ln ri7 
U.^DztU.U i 


i nn gqj-h m 
iuy .yoztu.ui 


195± 


81 


n m on 


n on 

u.zy 


U.oOztU.U f 


1 fl7 1 QJ.fl C\A 

1U f . lo±U.U4 


378± 


62 


T 1 n A O A 

HU4y4 


0.26 


U.2z±U.U0 


i nc coin in 


213± 


77 


T1 n^fic 
11UOOO 


n c.c 
U.oo 


Z.UOltU.lO 


iiri u 4_n nn 
11U.04±U.UU 


291± 


33 


T1 1 no^ 
lliuyo 


n Q7 
U.o / 


n Qo_i_n 1 n 
U.ozztU.lU 


i c\a OK4-H no 

lU4.ZDltU.UZ 


288±138 


T1 01 1 o 
llzllz 


U.OO 


i /i q j_n 1 fi 
1.4o±U.lo 


1 f\1 CfiJ-fl fiQ 
1U ( .OU±U.Uv5 


239± 


37 


112540 


0.60 


o c\o J_n 1 o 
2.Uo±U.lz 


i 1 n f;nj_n no 
11U.0U±U.U2 


243± 


21 


T1 QCQO 

iioooy 


n on 


U.o / ±U.Ud 


i c\a m J_n nc. 
1U4.U1±U.U0 


223± 


53 


T1 A Q/l fi 
114040 


n 

U.oO 


1 nfiJ-H on 
l.UOztU.ZU 


i nfi c; c_i_n nn 
lUD.OOztU.uy 


281± 


82 


T1 A Qfl/1 

114oy4 


n Q7 
U.o / 


n a e J_n i n 
U.4o±U.lU 


1 C\A QOJ.fl H7 

lU4.oo±U.U / 


325±113 


T1 K~\ Qfl 
llOloU 


n 7*3 
U. i o 


CA 

•CU.04 








T1 c,oc,n 

110ZOU 


n q.q 
u.oy 


n Qn4-n ns 
U.oUztU.Uo 


i nn o7J.n m 
iuy .z ( ztu.ui 


238± 


96 


T1 C/|fiO 
1104OZ 


n a o 
U.4z 


n oej_n r\Q 
U.^OztU.Uo 


i c\a QQj.fi nc 
lU4.5v5ztU.UO 


138± 


59 


116487 


0.27 


< 0.21 








117028 


0.40 


0.34±0.10 


104.30±0.04 


281±145 


117132 


0.52 


1.07±0.10 


109.67±0.04 


378± 


72 


117208 


0.69 


2.59±0.16 


110.54±0.02 


294± 


28 


118470 


0.37 


0.32±0.09 


106.86±0.01 


243±105 


119297 


0.27 


0.47±0.09 


106.20±0.11 


308± 


96 


120414 


0.29 


0.22±0.08 


106.10±0.11 


201± 


96 


122491 


0.33 


0.26±0.06 


107.03±0.02 


184± 


62 


123327 


0.33 


0.28±0.08 


104.06±0.09 


220± 


94 


123365 


0.42 


0.77±0.10 


108.30±0.03 


251± 


49 



arc the most luminous CO emitters in the local uni- 
verse, with their CO luminosity matching or exceeding 
those of the most luminous spirals. Unlike most late 
type galaxies that follow a nearly linear Lpjf{-L' co cor- 
relation (L FIR /Mn 2 =l-10 L & /Mq) however, ULIRGs 
are over-luminous in the far-/i? for their CO luminos- 
ity (Lfir/Mu 2 =10-100 Lq/Mq), leading to a significant 
up-turn in the Lfir-L' co correlation. 

Our new survey of 29 ULIRGs has increased the num- 
ber of ULIRGs with CO measurements to a total of 56. 
The resulting improvement in the total sample size and 
statistics provides a better definition of the L' co distribu- 
tion as shown in Figure 8. By focusing on the CO-bright 
and more luminous systems at slightly higher redshifts 
(z < 0.3), Solomon et al. (1997) established the devia- 
tion of ULIRGs from the late type galaxies in the field 
in the observed Lfir-L' co correlation. Our study has 
imporved the statistics particularly on the lower CO lu- 
minosity end of the nearby sample at z < 0.1. 

PG QSOs as a group have Lfir/L' co = 67±35 
L Q /[K km s -1 ], intermediate between ULIRGs (145±79) 
and late type field galaxies (22±21). As the result of 
probing the low CO luminosity of the local ULIRG pop- 
ulation, ULIRGs and PG QSOs now appear to con- 
nect more smoothly in the LpjR~L' co correlation plot 
(Fig. 8). The well known correlation between super- 
massive black holes (SMBHs) and their bulge stellar 
mass/velocity dispersion (Tremaine et al. 2002; Marconi 
& Hunt 2003) suggests that the formation/growth of 
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Velocity (km/ s ) Velocity (km/ S ) Fre ^ enc > 

Fig. 6. — Left) Comparisons of the 14m (dashed line) with 30m spectra (solid; Solomon et al. 1997). In the upper row, one example of 
CO bright (112540) ULIRG and one with intermediate CO line strength (100057) arc shown. In the bottom, each example show a CO faint 
object (110494) and one marginal detection (108572). The spectra from IRAM 30m is shown in mK in T m (, scale as presented in Solomon 
et al. (1997); Right) The RSR full-band spectra of the four ULIRGs are presented. The y-axis (in mK in TjJ scale) has been scaled down 
with the peak CO flux density of each case. The frequency ranges blocked with dotted lines are the regions compared with the study of 
Solomon et al. (1997) on the left-side. The upper arrow represents the optical redshift as Figure 3. 
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Fig. 7. — Distributions of the sample with different ranges of CO 
luminosities. The thick line represents 37 low-z ULIRGs (z < 0.3) 
from Solomon et al. (1997) and the shaded area represents 29 local 
ULIRGs (z < 0.1) from this work. 



SMBHs and their stellar hosts may be coupled. Pre- 
viously, some doubt was expressed on the ULIRG-QSO 
connection owing to the apparent discrepancy in their 
LfirI L'qq ratios (e.g., Yun et al. 2004). However, with 
a better agreement revealed by our new ULIRG CO sur- 
vey, this apparent discrepancy (driven by poor statistics 
and bias in the earlier surveys) is no longer a serious chal- 
lenge for the ULIRG-to-QSO evolutionary scenario (see 
5.2 for further discussion). 

A small number of ULIRGs with extreme ratios of 
Lfir/L' co > 250 Lq/[K km s _1 ] have been identified by 



our new ULIRG CO survey (Fig 8 and 9). PDS 456, the 
most luminous QSO in the local universe at z — 0.184, 
and an ULIRG at z = 0.059, 108572 were the only lu- 
minous infrared sources previously identified with such 
an extreme Lfir/L' co ratio (Yun et al. 2004). Yun et 
al. (2004) have speculated that these sources may be ex- 
tremely rare ULIRG-to-QSO transition objects based on 
the warm IR color and small Lfir/L' co ratio. Our 
study has identified seven additional sources, and they 
are discussed further below, in § 5.4. 

5.2. AGN Activities as the Origin of Large FIR-to-CO 
Luminosity Ratios for ULLRGs 

A number of LR luminous galaxies display characteris- 
tics of Seyfert activity in their optical spectrum, indicat- 
ing that dust heating by AGN activities may make an im- 
portant contribution to the LR emission. Although the 
actual frequency is poorly determined (Veilleux 2006), 
about 5% of LR bright galaxies with Lib, = lO lo ~ n £0 
and up to 50% of the sources with log Lm > 10 12 -^Lq 
show broad, high excitation lines characteristic of AGN 
activities (Veilleux et al. 1995; Kim & Sanders 1998; 
Veilleux et al. 1999; Kewley et al. 2001). An optical 
spectrum characteristic of Seyfert activity is found in 8 
out of 29 galaxies in our sample and 12 out of 29 in the 
comparison sample of local ULIRGs (Solomon et al. 1997; 
Gao & Solomon 2004a). Including those with a LINER 
spectrum (11 and 8), respectively) would increase the 
fraction of ULIRGs with AGNs. 

Given that the bulk of ULIRG luminosity emerges as 
dust-processed IR emission, optically thin tracers are re- 
quired to probe the presence of AGN activities with a 
greater accuracy. Also, since all massive galaxies are ex- 
pected to host a SMBH (e.g., Tremaine et al. 2002; Mar- 
coni & Hunt 2003), demonstrating that AGN activity 
is the primary source of luminosity is far more impor- 
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Fig. 8. — Lpm vs. L'qq of low-z galaxies with a range of far infrared luminosity. Spirals and starburst systems from Gao &; Solomon 
(2004b) and Solomon et al. (1997) are shown red triangles, z < 0.3 Paloma- (PG) QSOs from Alloin et al. (1992), Evans et al. (2001), and 
Scoville et al. (2003) are shown by green squares. ULIRGs by Solomon et al. (1997) are shown in blue circles. The yellow dots with a 
blue circle around represent the ULIRGs in our sample observed by RSR on the FCRAO 14m telescope. Non-detections are indicated with 
arrows. Dotted and dashed lines indicate constant ratios of Lpj _r/Mh 2 = 1, 10, and 100 using the standard CO-to-H2 conversion (Mh 2 
[Mq]= 4.6L' co [K km s _1 pc 2 ]). 



tant than establishing the presence of a hidden AGN. 
For example, Nagar et al. (2003) have reported that a 
high fraction (-75%) of ULIRGs with Seyfert or LINER 
emission lines host compact, non-thermal radio sources 
among the IRAS 1 Jy ULIRG sample. In contrast, the 
analysis of the Infrared Space Observatory (ISO) mid- 
IR spectra (Genzel et al. 1998; Rigopoulou et al. 1999) 
and IR photometry (e.g. Farrah et al. 2003) have shown 
that starbursts dominate the luminosity output in most 
ULIRGs while only 20-30% of ULIRGs are mainly pow- 
ered by a central AGN in the low-z Universe. 

Since the Lpm/ L' co may be an indicator of the ef- 
ficiency of converting gas mass into luminosity, we ex- 
plore whether this ratio may offer a useful insight on the 



nature of the powering sources for the large IR lumi- 
nosity of ULIRGs. In Figure 10, we plot mid-/i? color 
/25//6O1 as a proxy of luminous obscured AGN (de Grijp 
et al. 1985), against the Lpj^/ L' co ratio for the local 
ULIRGs, PG QSOs, and spirals. These three populations 
are fairly well separated in this figure with only a small 
overlap among them. Spirals and ULIRGs have similar 
mid-/i? colors, but ULIRGs are distinguished by their 
larger Lpm/ L' co ratios. PG QSOs show much warmer 
mid-/i? colors (/25//6O >0.2), and their Lfir/L' co ra- 
tios are comparable to or somewhat smaller than those of 
ULIRGs. Only -10% of the ULIRGs show mid- IR colors 
as warm as QSOs. The ULIRGs with Seyfert or LINER 
spectra are not clearly distinguishable by their mid-/i? 
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Infrared and CO Properties of the Sample 
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a The redshift measured by using the CO line. For non-CO detec- 
tions, an optically measured z is quoted and indicated with aster- 
isk. 
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Fig. 9. — LfirI ^co vs ' ^cc The same colored symbols are 
used as Figure 8. The dashed lines and the dotted lines represent 
L' CQ ±la (8.5±4.6x 10 9 K km s" 1 pc 2 ) and L FIR /L' CQ (145±79 
Lq [K km s _1 pc 2 ] -1 ) of 56 ULIRGs collected in this study. The 
error bars of our sample reflect the lc uncertainties in our CO 
measurements. 
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Fig. 10. — Infrared color (/25//60) vs. Lpjr/ L'qq. The same 
colored symbols are used as as Figure 8. The dashed line represents 
/25//6O = 0.2, the borderline of warm and cool ULIRGs. 



color. ULIRGs as a group are clearly distinct from PG 
QSOs in both /25//6O color and in their Lfir/ L'qq ra- 
tios. These trends imply that the SEDs of most ULIRGs 
are distinct from those of QSOs and that their mid-/i? 
emission may be indicative of dust heating by a starburst 
rather than AGN activities. Compared with non-ULIRG 
starburst systems, however, ULIRGs show a higher effi- 
ciency of converting molecular gas into stars as indicated 
by their high Lpm/ L' co as other studies have previously 
found. 

Another commonly used optically thin indicator of 
AGN activity is a compact, non-thermal nuclear radio 
source that produces radio emission in excess of the ex- 
pected trend from the well-known radio-/!? correlation 
("radio-excess" with q <C 2.1; see Condon 1992; Yun et 
al. 2001). For example, the majority of the high redshift 
QSOs detected both in mm/submm (rest frame far-/i?) 
show radio continuum well in excess of their iax-IR lu- 
minosity (Yun et al. 2000; Petric et al. 2003; Wang et 
al. 2007). A quantitative measure of the logarithmic 
FIR-to-radio flux density ratio, which is known as "g"- 
parameter, is derived using the following definition, 



log 



FIR 



3.75 x 10 12 W m- 2 



- log 



Si. 4 



W in" 2 Hz" 1 



(5) 



where S1.4 is 1.4 GHz radio flux density in W m -2 Hz -1 
(see Condon 1992). Lower q values correspond to higher 
AGN activity. In Figure 11 we plot g-parameter against 
Lfir/ L' co . As mentioned earlier, Nagar et al. (2003) 
have reported a high frequency (~75%) of compact, non- 
thermal radio sources among the IRAS 1 Jy ULIRG sam- 
ple with Seyfert or LINER emission lines. Among the 29 
ULIRGs we studied, however, only one object (114394) 
shows a clear radio-excess, and the radio-excess ULIRGs 
are rare among the local ULIRG population in general 
(see Fig. 11 and Table 4). In the local Universe, only 
about 1% of the 7i?-selected galaxies show clear radio- 
excess (Yun et al. 2001), and only about 10% of optically 
selected QSOs are radio- loud (Kellermann et al. 1989). 
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Fig. 11. — A plot of radio-to-FIR luminosity ratio ( "g-factor" ) . 
The same colored symbols are adopted as Figure 8. The long- 
dashed and dotted lines represent the mean q and the scatter of 
Condon (1992)'s sample of spiral galaxies (2.3±0.2). 
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Fig. 12. — L FIR /L' co vs. FIR luminosity. The dashed lines 
represent the typical star formation effeiciency (SFE) of the ULIRG 
sample (145 L Q /Li » 116 L©/M ) and L' co (8.47 x 10 9 L t ), 
respectively. Those of the spiral sample (22 Lq /L; 4-5 Lq /Mq , 
and 3.15 X 10 9 L;) are indicated with dotted lines (the conversion 
factors for different populations are given in § 5.4). 



Therefore, we can conclude that there is little evidence 
for a high level of radio AGN activity among the ULIRG 
population or alternatively can infer a low duty cycle of 
the radio AGN phase. 

5.3. Merger-induced Starburst Activities as the Origin 
of Large FIR-to-CO Luminosity Ratios for 
ULIRGs 

Since there is little direct evidence linking a large FIR- 
to-CO luminosity ratio and AGN activity, we next ex- 
plore whether the merger-driven starburst phenomenon 



can account for the large luminosity ratios observed. In 
Figure 12, we plot the Lfir/ L' co ratios of ULIRGs and 
non-ULIRG systems as a function of Lfir- A typical 
star formation rate of nearby CO-rich spirals (4 LqMq 1 ) 
is similar to that of galactic clouds and is generally at- 
tributed to recent formation of OB stars in a quiescent 
disk (e.g. Solomon & Sage 1988). The Lfir/L' co ratio 
of 20 LqMq 1 , shown by a dotted line, is more typical 
of local starburst galaxies. Yun et al. (2004) and others 
have noted that PG QSOs display a similar Lfir/ 'L'co 
ratio and have proposed a closer physical and evolution- 
ary link with the local starbursts, rather than with the 
ULIRGs. Most local ULIRGs have a Lfir/ L' co ratio 
well above the value of non-ULIRG starburst systems, 
implying that the ULIRG phenomenon not only requires 
the largest amount of molecular gas found in the most 
gas-rich systems in the local universe, but it also requires 
a much more efficient mechanism that converts the gas 
mass to luminosity. Another potentially important clue 
to the underlying physical mechanism is that the disper- 
sion in the F I R-to-CO luminosity ratio for the ULIRGs 
is the largest among the different populations compared. 
This dispersion may reflect varying evolutionary stages 
of the ULIRG phase. 

A useful insight on the observed Lfir/ L' co ratio (also 
dubbed "star formation efficiency" Young & Scoville 
1991) for the starbursts and ULIRGs is offered by nu- 
merical studies of merging disk galaxies. Motivated by a 
high frequency of morphological peculiarities associated 
with ULIRGs (Sanders et al. 1988a; Farrah et al. 2001), 
Mihos & Hernquist (1996) have investigated the response 
of interstellar medium (ISM) during a merger of two gas- 
rich spirals and the subsequent secular evolution of the 
triggered starburst activity. A particularly relevant re- 
sult from this study is that the compression and concen- 
tration of gas (details of which are shown to be highly 
dependent on the initial conditions) during the course of 
the merger lead to spikes of elevated starburst activity 
lasting 5-20 xlO 6 yr (see their Fig. 5). In a more recent 
simulational study, Matteo et al. (2008) have investigated 
the enhancement in star formation during galaxy inter- 
actions, suggesting that major mergers can increase the 
star formation efficiency (SFE) by a factor of a few tens 
compare to that in unperturbed disks which appears to 
be consistent with what is shown in Figure 12. 

We interpret this prediction using our data by drawing 
lines of constant CO luminosity L' co (long-dashed lines) 
in Figure 12. The long dashed line on the right side cor- 
responds to the average CO luminosity of all ULIRGs, 
and it goes through the upper envelope of the CO lumi- 
nosity associated with CO-rich spirals in the field. While 
the long dashed line on the left corresponds to a L' co 5 
times smaller, approximately the L* value of the local 
CO luminosity function (10 9 K km s— 1 pc 2 ; Keres et al. 
2003). 

As one of the CO-rich field galaxies experience a large 
spike in the massive star formation rate (induced by a 
galaxy-galaxy collision in the cases of the Mihos & Hern- 
quist study), both Lfir and Lfir/ L' co ratio would in- 
crease along the right long-dashed line, reaching the area 
occupied by the ULIRGs. As the starburst activity fades, 
the galaxy would climb back down along the same line. 
Therefore, the apparent broad trend between Lfir/ L'co 
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TABLE 4 

AGN Diagnostics for the Extreme L fir /L' co ULIRGs 
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cf The mean ratio of L FIR /L' CQ of the ULIRG, the QSO, and 
the spiral sample collected in this study are 145 (113 excluding 9 
extremers in the table), 67, and 22 Lq/Lj, respectively. 

ratio and Lfir seen in Figure 12 can be naturally ex- 
plained by secular changes in the starburst luminosity 
(externally triggered or internally by dynamical instabil- 
ities such as bars or spiral density waves) to the gas mass 
ratio. 

One may expect the evolutionary trajectory to be 
somewhat skewed from the right long-dashed line since 
the starburst activity would consume some of the gas and 
would possibly remove additional gas through energetic 
feedback processes (e.g., Heckman et al. 1990; Martin 
2005). On the other hand, the resulting higher gas pres- 
sure will lead to an increase in CO excitation and line 
luminosity, partly compensating for the gas mass loss — 
see the discussions of elevated CO luminosity in nuclear 
starburst regions by Solomon et al. (1997) and Downcs 
& Solomon (1998). These competing effects should con- 
tribute to the tightness in the observed correlation. 

5.4. Extreme Lfir/L' co Objects and ULIRG Evolution 

An important outcome of our ULIRG CO study is the 
discovery of a special group of ULIRGs; 110494, 105189, 
122491, 108572, 120414, 110173, and 115250 (in the or- 
der of increasing L' co ). These systems show extreme 
LfirI L'qq ratios which ranges from 250 to 310 Lq/Li 
(Fig. 8), which is greater than the mean of the rest 
ULIRGs by a factor of two or more. With the standard 
conversion factor for spiral galaxies, this corresponds to 
a star formation efficiency of > 50 Lq/Mq as shown 
in Figure 8, and potentially > 300 Lq/Mq, adopting 
the conversion factor for starburst systems (Downes & 
Solomon 1998). All but one (115250) had been previ- 
ously observed in CO by others (Table 1) and five were 
detected but 120414 although they have not been rec- 
ognized as a distinct group before. The only other ob- 
ject previously known for such a large Lfir/L' co ratio 
is PDS 456, the most luminous QSO in the local uni- 
verse (z < 0.3 Yun et al. 2004). These objects are read- 
ily identifiable in Figures 8, 9, and 12 by their extreme 
LfirI L'qq ratios, apart from the rest of the ULIRGs. 
Along with PDS 456, these seven ULIRGs have some of 
the lowest far-/i? luminosity, suggesting that a compar- 
atively larger fraction of their total IR luminosity arises 
in the mid-/i? wavelengths for a given Lfir/ L' co . While 
the extreme Lfir/ L' co of these systems is primarily at- 
tributed to low L' co , there are also two objects with a 



similar Lfir/ L' co but higher L'qq among the sample of 
Solomon et al. (1997), 103158 and 14070. These are two 
of the highest luminosity ULIRGs known, with nearly 10 
times larger far-/i? luminosity than the seven extreme 
ULIRGs identified in this study such that these may be 
a different class of objects altogether. However, this lu- 
minosity distinction may simply arise from the redshift 
restriction imposed by the RSR spectrometer. Our sam- 
ple includes the lowest redshift ULIRGS corresponding 
to lower luminosities. 

As a first step toward understanding the mechanism(s) 
responsible for their extreme Lfir/ L' co ratio, we ex- 
amine in detail the known properties of the individual 
sources, including the two from the Solomon et al. (1997) 
sample: 

• IRAS 03158+4227. This z = 0.134 ULIRG shows 
a second core which is much fainter than the one 
in the center. While the optical morphology is sug- 
gestive of a merging event, this object is found near 
AGN dominated systems in Spoon et al. (2007) 's 
diagnostic diagram for mid-infrared spectra of IR 
galaxies, with a relatively small equivalent width of 
the 6.2 /iin PAH emission feature and a weak 9.7 
/j,m silicate strength. 

• IRAS 14070+ 0525. This compact z = 0.264 
ULIRG does not show any obvious signatures of 
interactions or merging. This system however, has 
disturbed isophotes (Veilleux et al. 2002) and mor- 
phologically classified as a remnant by Dasyra et 
al. (2006b). Its spectrum shows characteristics of 
Sy2 (Veilleux et al. 1999). 

• IRAS 05189-2524. The optical image of this 
z = 0.043 ULIRG suggests that it is in a late stage 
of a merger with a highly obscured, compact single 
nucleus (Farrah et al. 2003). It is one of the six 
ULIRGs in the IRAS Bright Galaxy Sample with 
the warmest /25//6O color (Sanders et al. 1988a). 
Presence of both an obscured AGN (Young et al. 
1996) and a buried starburst (Dudley 1999) is sug- 
gested by a Sy2 spectrum (Veilleux et al. 1995) 
and an X-ray spectrum characteristic of a Comp- 
ton thin AGN and a thermal component (Risaliti 
et al. 2000; Lutz et al. 2004). By modeling the 
observed SED, Farrah et al. (2003) conclude that 
> 30% of the total IR luminosity is contributed by 
the AGN. 

• IRAS 08572+3915. This z = 0.059 ULIRG shows 
two cores which are well separated by 6 kpc. Both 
the stellar nuclei are dominated by Ha emission 
(Colina, Arribas & Monreal-Ibero 2005). It is also 
a well-known warm object (/25//6O = 0.23) with 
a strong evidence for an AGN (Soifer et al. 2000; 
Farrah et al. 2003). 

• IRAS 10173+ 0828. This z = 0.049 ULIRG looks 
disky compared to most ULIRGs with a highly 
asymmetric outer disk. It is a OH megamaser 
galaxy (OHM), which is thought to be caused by 
infrared radiation from the surrounding environ- 
ment (Baan 1985). Most OHMs are known to 
be warm (U)LIRGs (Chen, Shan, & Gao 2007) 
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although this system is not a particularly warm 
ULIRG among the sample with /25//6O ~ 0.11. 

• IRAS IO494+4424. This z = 0.092 ULIRG 
presents two spikes to the north-northwest which 
(Veilleux et al. 2002) suggest tidal origin. In 
their SPITZER spectroscopic study, Imanishi et al. 
(2007) find no signature of a buried AGN in this 
system, which is consistent with the SED study by 
Farrah et al. (2003). The upper limit of the AGN 
contribution to the total IR flux estimated by Far- 
rah et al. (2003) is less than ~L3%. 

• IRAS 15250+3609. This z = 0.055 ULIRG shows 
a ring-like structure around a spheroid in the cen- 
ter (Scoville et al. 2000) with a much dimmer nu- 
cleus 0.7 arcsec away from the center (Farrah et 
al. 2003). Lutz et al. (1999) have classified it as 
a starburst/LINER based on its mid-IR spectrum. 
However, in Farrah et al. (2003) 's SED study, this 
object has been found with one of the largest AGN 
contribution (Lagn/ Lir si 0.6) among the sample. 

• IRAS 20414- 1651. This z = 0.087 ULIRG was 
observed in CO by Mirabel et al. (1990) but not 
detected. It shows a main nucleus with highly elon- 
gated isophotes (Veilleux et al. 2002). Murphy et 
al. (1996) find a double core in the optical which 
appears to be one system in the infrared. Farrah 
et al. (2003) find its AGN contribution to the total 
IR luminosity of ~ 13%, below the mean of their 
sample (23%). However, Dasyra et al. (2006b) find 
a black hole mass of > 10 8 Mq (larger than that 
of 105189 by a factor of w 3.5) based on the bulge 
dispersion relation (Tremaine et al. 2002). 

• IRAS 22491- 1808. This z = 0.078 ULIRG is a 
close (2.2 kpc) pair with two tidal tails (Veilleux et 
al. 2002). Cui et al. (2001) have proposed a multi- 
ple merger origin for this system, which is also sup- 
ported by its mid-/i? spectrum (Lutz et al. 1999). 
However, this ULIRG is one of the few objects with 
the largest AGN contribution (Lagn/Lir ^ 0.7) 
in the sample of Farrah et al. (2003) 's SED study 
along with 115250. 

Based on their high Lfir/L' co and evidence for 
AGNs (except 103158 which has not been well studied, 
and 110494 with the lowest Lfir/ L' co ), we speculate 
that these ULIRGs arc currently consuming enormous 
amounts of molecular gas to feed the central black holes, 
and are likely to become QSOs. The rarity of PDS 456- 
likc QSO however, suggests that the timescale of forming 
classical quasars in this way is relatively short compared 
to the timescale of growth of an AGN in a ULIRG phase. 

The interconnection between ULIRGs and QSOs has 
been suggested by Sanders et al. (1988a, b) who have 
demonstrated the similarity between the SEDs of warm 
ULIRGs (/25//6O > 0.2) and those of QSOs. Warm 
ULIRGs also morphologically resemble QSOs with a 
more prominent spheroid, weaker tidal features, and 
brighter nuclei compared to their cooler counterparts 
(Veilleux 2006). These observations suggest that cool, 
starbursting ULIRGs may go through an AGN-like warm 
ULIRG phase and eventually become optically selected 



QSOs once the burst of star formation decays and the 
nucleus sheds its obscuring dust (Sanders et al. 1988a; 
Yun et al. 2004; Veilleux 2006: Dasyra et al. 2006a). 

The transition from a ULIRG phase to a QSO phase 
can be driven by rapid conversion of gas into stars and/or 
a subsequent growth of a supermassive black hole (Yun 
et al. 2004). Due to rapid consumption of gas, transi- 
tion objects are expected to deviate from the standard 
ULIRG L fir-L'co correlation in the way that results in 
a larger L fir for a given molecular gas mass compared to 
most ULIRGs. There are only a few such ULIRGs that 
have been found in previous CO studies which are bi- 
ased to luminous ULIRGs and one QSO, PDS 456 (Yun 
et al. 2004). In our study however, we find four addi- 
tional ULIRGs with PDS 456- like L FIR /L' co along with 
108572 which already has been known (Gao & Solomon 
1999). 

6. SUMMARY 

We have presented 12 CO J = 1 — > observations of 
29 ULIRGs at z=0.043-0.11, using the Redshift Search 
Receiver (RSR). In total 27 systems have been de- 
tected including 9 new detections, which has increased 
the number of local ULIRGs with CO measurement 
by 15%. The CO line luminosity L' co of our sample 
ranges from 1.2 to 15.3 x 10 9 K km s _1 pc 2 . Adopt- 
ing the CO-to-H2 conversion factor for ULIRGs/QSOs 
(M H2 =0.8 L' co M q [K km s" 1 pc 2 ]" 1 ), the inferred cold 
gas mass in H2 form in those ULIRGs is 1 - 12xl0 9 Mq. 
We have investigated L' co of local ULIRGs as function 
of IR and radio properties. 

1. ULIRGs are 10-100 times higher in Lfir/ L' co 
than spirals/starburst galaxies and form a contin- 
uous track with low-z QSOs, consistent with re- 
sults from previous studies. However as a result 
of probing deeper in CO faint objects, our survey 
finds a broader range of Lfir/L' co , resulting in a 
smoother transition between the ULIRG and QSO 
populations. 

2. ULIRGs are well separated in /25//6O vs. 
LfirI L'co s P ace from spirals/starburst systems by 
much larger Lfir/ L' co . They are also distinct 
from QSOs by their cooler mid-/i? color (/ 2 5 /fan < 
0.2). 

3. The radio excess, q of ULIRGs is not significantly 
different from those of non-ULIRG populations. 

4. Seven ULIRGS with extreme values of 
Lfir/ L' co ^ 250 are identified. These are 
similar to the luminous quasar PDS 456. 

From these, we conclude that the power sources of 
most local ULIRGs are mainly sporadic starbursts which 
are likely to be driven by merging events, rather than 
AGNs. However, we always find some evidence for a 
powerful central source in those ULIRGs with extreme 
Lfir/ L' COl which may be representative of transition 
objects to QSOs. These objects are not distinct from 
most ULIRGs in other properties which might be due to 
a short timescale of the transition. 
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